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Abstract 



The 0(N) non-hnear sigma model in a L)-dimensional space of the form ^ x T*"^, 

^ \ R-^^^^ X S*^, or T^^ X is studied, where R*^, T*^ and S*^ correspond to flat space, 

' a torus and a sphere, respectively. Using zeta regularization and the expansion, the 

\ corresponding partition functions — for deriving the free energy — and the gap equations 

■ are obtained. In particular, the free energy at the criticial point on R^'^+^ x S^^"^^ vanishes 

' in accordance with the conformal equivalence to the flat space R'^. Numerical solutions 

Q^. of the gap equations at the critical coupling constants are given, for several values of D. 



The properties of the partition function and its asymptotic behaviour for large D are 
discussed. In a similar way, a higher-derivative non-linear sigma model is investigated 



^ . too. The physical relevance of our results is discussed. 
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I Introduction 



Zeta regularization |l| (for a review see 0) is a very powerful and elegant method for reg- 
ularizing the divergences that appear in quantum field theory. It has found lots of different 
applications, from the calculation of the vacuum energy density or Casimir energy correspond- 
ing to very different configurations (different fields, different spacetimes, different boundaries) 
to its application in wetting/ nonwetting phenomena in actual condensed matter and solid state 
systems, to the analysis of phase transitions coming from the study of effective potentials in dif- 
ferent context, topological mass generation, Bose-Einstein condensation phenomena, evaluation 
of the partition function in string and p-brane theories, etc. f^. From a more mathematical 
point of view, the method has allowed the computation of the basic operator tr log (□ +X) on 
a curved manifold, that is very important in quantum gravity and cosmology. 

In this paper we will study the zeta-function regularization of the 0{N) non-linear sigma 
model in an arbitrary number, D, of dimensions and on spaces of the form R^^^^ x T^^, R^^*^ x 
S^^, and T*^ x S^. In flat spacetime this model has a big number of different applications (see 
for a review), in particular, in the theory of critical phenomena and solid state physics. For 
example, such a three-dimensional model, which is known to be renormalizable in this case, 
may be used in condensed matter physics as an effective field theory of the two-dimensional 
quantum antiferromagnet [O. In particular, the 0{N) non-linear sigma model on x R^ 
may be applied to describe the low-temperature properties of the quantum antiferromagnet 
T8| . Recently, in connection with the study of higher-dimensional conformal theories 0] — 



where one can use well-known 2d conformal field theory techniques — such a model has been 
considered in three dimensional curved spacetime [^. The critical properties of the model were 
also discussed. 

Here we will extend this important analysis, by studying the 0{N) non-linear sigma model 
in topologically non-trivial spaces of constant curvature in arbitrary dimension D. Using zeta- 
function regularization and the expansion techniques we will obtain the partition function 
and the gap equation in each of the cases considered, and also numerical solutions of the gap 
equations for some spaces. The asymptotic behaviour in the limit D ^ oo will be considered 
too. Finally, a higher-derivative generalization of the 0{N) non- linear sigma model will be 
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introduced and the corresponding zeta-function on a flat but topologically non-trivial space 
will be obtained. Numerical solutions of the gap equations in this last model, for different 
values of will be discussed too. The relevant partition function and gap equations are 
presented in sects. 2 and 3, while the applications to R^"*^ x T^^ and R^^*^ x S^^ appear 
in sects. 4 and 5, respectively. Our higher-derivative model is introduced in sect. 6. In the 
conclusion we comment on the relevance of our results and discuss future perspectives. The 
appendix is devoted to the calculation of the necessary zeta functions for the spaces under 
consideration. 

II The 0{N) nonlinear sigma model in D dimensions: 
partition function 

Let us consider an arbitrary D-dimensional space of constant (or zero) curvature. It is well- 
known that the scalar conformally invariant D'Alembertian operator on such manifold is 

-a + ^R, (II.l) 

where ^ = 4^Eij- In what follows we use Euclidean space notations. The partition function of 
the 0{N) nonlinear sigma model in D dimensions we will be interested in is given as follows: 

Z[g] = I V^Vaexp J d^'x^ [4>\-n + ^R)^^ + '-!)]}, (n.2) 

where 0* are scalars in curved spacetime, i = 1, . . . , N, ^ = ^1^^^-^ is the conformal coupling, a 
is an auxiliary scalar introduced in order to keep the constraint = 1, coming from the 

condition of 0(A^)-invariance, and A is the coupling constant. Note that a has no dynamics, as 
it just plays the role of a Lagrange multiplier. Observe also that we have choosen to work with 
the conformally invariant operator — □ + ^i?, in order to better understand the conformally 
invariant properties of the model. 



The theory (|II.2|) was extensively studied in in three-dimensional curved spacetime at 
its nontrivial fixed point. In particular it was shown, in the large- limit, that such a model is 
an example of a conformal field theory at a nontrivial fixed point. Investigation of this theory 
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in different spaces of constant curvature has suggested that what distinguishes a given model 
is not curvature, but the conformal class of its metric. 

Our purpose here will be to study the theory in D-dimensional curved spacetimes of con- 
stant curvature near the nontrivial fixed point in the 1/A^-expansion. What is even more 
important, we will analyze, in addition to the large- iV limit, some situations involving the limit 
D —>■ (X), which is relevant to dimensional dependence investigations in a number of quantum 
systems [|], 0. Numerical solutions to the gap equations will also be given. 

It is convenient to rescale \fX(j). Then 



Z[9] 



V(j)Va exp 



where the mass- dimensions of the fields and parameters are 



D 



(J 



D-2, 



[e] = 0, 



(11.3) 



(11.4) 



and the dependence of Z on the metric (yf^,^ is explicitly shown. Note that sometimes it is 
convenient to rewrite the partition function ( [11. 3| ) as explicitly regularized. In particular, if one 
uses a cutoff A for regularizing, it may be adequate to do the change 



A(A) 



in order to 



work with a dimensionless A (A). 

The aim is to study the above theory in the large- limit keeping, as usual, A^A fixed 
as A^ oo. The spacetime dimension D will be arbitrary. Integrating out the first A^ — 1 
components of and rescaling the A^th component 0Ar to 

A^ - 1 

{Trlog(-n + ^i?+(T) 

(11.5) 



N-l 



■>N, 



and ^^^^A to A, we get 



z[g] 



l)(j)]\rT)a exp 
1 



+ - J d'^x^ 



a 



-□ + ei? + cr)0^- — 



III The gap equations 



Since we shall deal with manifolds of constant curvature, we will look for a uniform saddle 



point: a{x) = m?, (f)N{x) = b. Extremizing the action (|II.5| ) with respect to 0Ar(a;) maintaining 
a{x) fixed, and the other way around, we obtain the gap equations 

(-□ + ^i? + m2)6 



G{x, X] m^, g) + h'^ — ^ 



0, 
0, 



(III.l) 



where 

G{x,x;m'^,g) = {x\{-a + + m^)-^\x) (111.2) 

is the two-point Green function at equal points. Once the sohitions to these equations have 
been found, it is sensible to evaluate the free energy density W at the saddle point to leading 
order in 



W[g,\] 



N 
~2 



.m 

T 



d^x^- 

Applying zeta-function regularization one defines 

G{x,x]w? , g) = YimGix.x-.m^ , g] s) oc lim(a;|CA^(s)|a;), 

s^l s—*l 

Tr log(-n + e/? + m2) (X -C^wW- 

Cm{.s) is the spectral zeta function of the operator — □ + + rn? on the spacetime under 
consideration, i.e. 



(III.3) 



(111.4) 



Cm{s) 



Tr {-U+^R + m^ 



(III.5) 



The proportionality factors — not explicitly written in ( [1II.4| ) — are determined, in each case, 
by the normalization of the physical states, and they have such a form that the dimensionalities 
match. 

By heat-kernel series analysis, it is known that the short-distance divergences of this two- 
point Green function depend in general on the curvature of the spacetime, except for the leading 
pole which is independent of g and, therefore, present in all cases. As a result, in flat spacetimes 
with nontrivial topology such as or R^"^-'^ x T^^, this singular behaviour will be the same as 
in the flat space R^. That is why it is natural to study R^ first. According to this reasoning, 
the critical value of A — at which the theory becomes finite — will be the same for all these 
spacetimes. 

The gap equations in the spacetime R^, after momentum cutoff regularization, read 

m^h = 0, 
= 



A 



A(A) 



— Ga(x, x: m 



(111.6) 



Here GA(x,x;m'^,g) means the Green function obtained when setting a cutoff A on the norm 
of the integrated momentum. Studying their solutions, one finds that, for 6 = m = 0, 



A 



Ac(A) 



(A) d'^k 1 



A 



D-2 



(27r)^P ^4^)B/2r f - r 
5 



if D > 2, 



(III.7) 



and, for 6 = 0, m 7^ 0, 
A A 



Ae(A) A(A) 



(III.8) 



where e(A) are terms vanishing as A — oo. Since T [1 — < for odd D > 1, this indicates 
the unphysical character of the solution when A < Ac. A more detailed study shows that A = Ac 
is a critical value separating two different phases. 

In the same spacetime, when using zeta-function regularization the second gap equation 
becomes 



lim 



1 



where 



G{x, x; m^, g; s) 



X{s) 



G{x, x; m^, g; s) 



1 



^D-2s r (s - f ) 



(III.9) 



(111.10) 



(27r)^(P + m2)^ (47r)f r(s) 
Following the line of thinking of ref 0, we set the values m = 6 = from the discussion in 
cutoff regularization, and realize that now the only consistent way out is 

1 



lim 



0, 



>i Ac(s) 

which gives the critical value of A in this regularization. 

In curved spacetimes without boundaries, heat-kernel expansion gives 

/oo „ 



(III.ll) 



(47r)«/2 
1 



(47r)^/2 _ 



^ - 1 
L 2 



^ -2" 
2 ^ 



OD-I T2 
1 
2 



+ /i(A), 



for odd D, 



^0 ~xD-2 



D 



A' 



+ 



D 



«2 



-2 



+ . . . + aD-iA."^ + aD-2 log A^ 



+ /2(A), for even D, 

(III.12) 

where /i(A) and /2(A) are terms that become finite when A — 00 and where, for consistency, 

D-2 

we have taken A^ = F (^^^ ^ A^. As usual, the a2n's stand for the even Seeley-Gilkey 
coefficients ^ , and we have taken into account that a2„+i = in the absence of boundaries. 
Since oq = 1, for the first D's we immediately get 
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D 



3 
4 
5 
6 
7 



divergences of Ga{x, x-.rn?,. 



1 

(47r)3/2 
1 

(47r)2 
1 



2A 

+ a2 log A^ 
(|A3 + 2a2A 



(47r)5/2 

^ (|A3 + a2A2 + a4 logA^ 



(4 



(fA^ + fasA^ + 2a4A 



(47r)7/2 

The D = ?) result does not depend on 02, 04, . . . and is thus independent of the curvature. In 
consequence, it is the same for any 3-dimensional space without boundaries and therefore it is 
enough to find it in R^. For conformally flat manifolds of any D, the property i? = makes 
02,04, . . . vanish with the same consequence, i.e. , we can get by with the critical values in 
R^. This fact will be used in the study of R^"*^ x T^. 

When i?7^0inD>3,a particular study of every particular situation is called for. In 
spaces of constant positive curvature, such as spheres or some products of R"^ by spheres, the 
first gap equation reads {^R + m^)6 = with i? > 0, and admits no other solution than 6 = 0. 
From the second equation, = Ga(x, x; m^, (7), where the divergent parts of the r.h.s. as 

A ^ 00 are given by (pi.l2| ). Next, we form the difference between the m = and the m > 
cases: 



A 



A 



A^=o(A) A(A) (47r)^/2 



02(771 = 0) — 02(m)^^_4 ^ 04(m = 0) 



^ -2 
2 ^ 



^ -4 
2 ^ 



(II 



.13) 



To obtain 02„(m) from 02n(m = 0) is just a matter of replacing ^ — ^ + ^ in every expression 
of these coefficients. As an example we consider D = 5, whose divergences appear in the 
preceding table. Since 02(m = 0) = — |) R, 02(m = 0) — a2{m) = m? independently of the 



value of R. After dividing by A we get 

1 



3/4 



-2m^ 



(III.14) 

A. 



A„=o(A) A(A) (47r)5/2 

This equality stops making sense if A < \m=Q, which can be interpreted by regarding Am=o 
as the critical value separating two phases. Since this was obtained by setting 6 = and m = 
those are the values we shall set for computing ^ in zeta regularization [^. 
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IV R^-^ X 



Putting the adequate normalization factors, the Green function and the free energy at the 
critical value of A read 

T7 = ~o Ti" Cro-m^t^m (0). 



N 2 V27r. 

After calculating the zeta function for our operator on this spacetime (see the Appendix), we 



can write 



' / D-M pm\ /pm\^-2« / D 



(rv.2) 

The second term tells us that the singularities at s = 0, 1 arc present only when D is even, 
independently of M. Im (s — ^) , is an integral of the type 



(IV.3) 



_ \ 

where Q{x) is the Jacobi function 

Q{x)^ e"'''^"'. (IV.4) 



oo 

,2 



The key point is that the integrand is well-behaved around i = 0, causing no new pole at 
z — M/2 ( i.e. at s = D/2). Further, Im can be expanded into a Dirichlet series, which reads 



where K is the modified Bessel function and \n\i stands for the Euclidean norm of n = 
(ni, . . .,ni). 

The particular cases we will calculate are those of odd D = 2d + 3, M = 1. Under these 
conditions the Bessel functions involved are just 

{d + k)\ 1 



Kd 1 (x) — K,,i (x) — \ — e "^y^ , , , - , 
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and we may interchange the summations to obtain 



2d 



d+3/2 f^P^^ 



2d+l 



)d+2 



Cr2<*+2xTiC1 ) — ^ 

-J^ {d + k)\ Li,+i+fc(e"^'") 



+ T -d 



<^+l Mil \ U\\ T i ^^—pm\ 



2d+3 



(IV.6) 



2d+3 (^+^ + ^)' Lirf+2+fc(e 



-^rf-d-- 

■'^k\{d+l-k)\2'' {pmY+'^+^ V 2 



where Li is the polylogarithm function. Setting c? = = 3), these formulas reproduce — as 
should be expected — the ones in [Q. 

We look at the solutions of the gap equations for the critical value y- = 0. The chances 

are: 



1. m = 0, 6 7^ 0. This leads to 



d+3/2 



2d+3 .25TT4v^^C(2rf+l). 



When (i = the r.h.s. diverges and no solution can exist. For d > 0, C^{2d + 1) > and 
there is a sign conflict which prevents the appearance of any solution on this side. 



2. m > 0, 6 = 0. Now we are posed with solving 

[d^k)\ 



E 



'^k\{d-k)\2' 



:(pm)'^-^Li,+i+,(e-^"^) + (pm) 



2d+\ 



0. 



(IV.7) 



Only (i = admits a relatively simple analytic solution because Lii is the only polylogarithm 
which can be trivially expressed in terms of elementary functions. In that case one gets (see 
also H) {pm)c = 21ogr, with r = , i.e. {pm)c — 0.9624. Of course, it is also possible 
to find this same value by numerically solving the above equation for d = 0, and this is precisely 
what we do for the next d^s. Afterwards, the value found is replaced into ( |IV.6| ) and ( [IV. 1| ) , so 
as to find the free energy W. We thus arrive at: 
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d 


{pm)c 


^(pm = (pm)J 





0.9624 


-0.1530 


1 


no solution 




2 


2.1775 


-0.0441 


3 


no solution 




4 


3.5504 


-0.0561 


5 


no solution 




6 


3.6841 


-2.2634 



The figure for ^ when c? = coincides with the numerical value of — -^ClS), derived in [^] with 
the help of polylogarithm identities from [|13|. The rest of the values are the first (and possibly 
only) solutions found after scanning a reasonable positive range. 

We can now study the asymptotic behaviours of these expressions for M (and D) and/or 
a (= pm/(27r)) going to infinity. Two cases will be considered: (i) a ^ 1 and M bounded, 
with M — D fixed, and (ii) a ^ 1 and M ^ 1 with a/M ^ const., again with M — D finite. 

Let us start with the first case. As everywhere the dependence on M is through D — M, it 
is enough to study Im{t, a), where r = s — {D — M)/2 with s = or s = 1. From the behaviour 

1/2 



K 



D/2-S 



(27ra;|f2|) 



TT 



4:iTa\n\ 



-27rQ!|n| 



(IV.8) 



we easily obtain that 



M 

JM(r,a)~vr^'^/^«(^-^)/2-E 

1=1 




2' T.i^'^ (p- 



2TTa\ 



(IV.9) 



Where Lip^(x) denotes the generalized polylogarithm, for which we have 



Li^pVx) 



E 



X 



X 



ni,...,ni=l ^fif + ■ ■ ■ + nf 

Taking this into account, we obtain 

/ D-M 
hi s — ,a 



+ 



X 



/I+3 



.VTT3 



X -Cl. 



< 27r^-^/2^(^-^)/2-s^2g-2™_ 



(IV.IO) 



(IV.ll) 
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We can now consider M to be large, of course, but never competing with a or the above 
expression loses its sense. In order to deal with both limits at the same time, we must consider 
the case (ii). From the well-known behaviour of Kyiyz) for constant z and v — >■ oo, by calling 

'27T(y\/~l j xii 

Ui = lim — — = const., Vi = r]{ui) = Jl + uf + log . = const., (IV. 12) 

a,D^oc D/2 -s 1 + ^1 + ^2 

we get 

V R^-^^ X 

The Green function and the Tr log contribution to the free energy are now 

G{x,x;m^,g;l) = ^^^^ d^m^m Cr^-^^ x ( 1 ) , (V.l) 

and 

respectively, where a is the radius of the sphere. Here we will consider the massless case only. By 

our discussion on gap equations, in D = 5 — = G{x,x; 0,g, 1). Therefore, in these conditions, 

Ac 

Crs-^xs'*^ (1) for m = tells us the critical value of A. Another reason for the m = choice, 
apart from simplicity, is that the M = D — 1 case is conformally equivalent to R"^ — {0}, and 
then, m = corresponds also to solutions for critical A. 

Taking the conformal coupling and the known form of the Riemann scalar for the sphere, 
we follow a method analogous to ref. |jTl| (see also |T^) to construct the required zeta function, 
which is 

r (.^ o.,2.-D+A/ ^ ^(^-^) , f S.(P,2.-D + M), forM = 2p + 2 



r(s) I S/3(p,2s-D + M), forM = 2p+l, 

(V.3) 



where 



Sa(p,2z) = ^^^^ X](-l)'«^'b-l)C// (-2p + 2A; + 22;-l,i 
S/3(P,2^) = |:(-l)'/5fc(p-2)Cij(-2p + 2A; + 2^), 



(V.4) 



(2p) 



fc=0 
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C,H and C,R denote the Hurwitz and Riemann zeta functions, and the ak and j3k coefficients are 

(V.5) 



«o(i) = 1, afe(i) = ^ + ^) ■ ■ ■ (^^ + ^) ' ^ - ^' 

/3o(j) = 1, = E + + A;>1. 

0<ii<...<ife<j 



Our notation is just shghtly different from that in ref. [|ri|]. In fact a^ip — 1) = -^a^'^{p — 1) 
and (3k{j) — 2) = — 1), where CC stands for the coefficients employed in [|TT|. These 

coefficients enable one to write 



fcmax{M) ^ M — 1^ M-1-2A: 



d{M,l)= E {-lfAk{M)[l + ^^) , (V.6) 

where = 2/{2p + - 1) and k^^^{M) = p for M = 2p + 2, while Ak{M) = 

2/{2p)\Pk{p - 2) and k^^M) = p - 1 for M = 2p + 1. 

In order to study the Green function and the Tr log contribution to the free energy we 
need to evaluate the above zeta function at s = 1 and its derivative at s = 0. Such quantities 
will be expressed in terms of 

^a{p,2z) = ^(_i)^«,(p-l)G(-2p + 2A; + 2z-l,i), 



^ (V.7) 



V(P,2z) = tAt E(-lMb-2)C;i(-2p + 2/c + 2 
{■^Pj- k=0 

and the results will follow Then, we have the following cases 
1. Even D-M = 2q 



(a) M = 2p + 2 

(g- 1)! 



i + 7 + ip{q) + 2 log a/i 



S,(p,-2g + 2) 



+2J:'^{p, -2q + 2)} + 0{e) 
+ ip{q+l) + 
+2E'Jp,-2g)}. 



Cr2,xs2p+2'(0) = 7r«^il!2a-2« {[7 + + 2 log a/i] S„(p,-2g) 



(V.8) 

As usual, /i is a parameter with mass dimension, introduced by redefining Cj^^^^{s) 
as /i^^^Cl^^^''^''^'^ (s), which renders the log arguments dimensionless. 
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(b) M = 2p + 1 

Cr2,xs2p+i(1) = vr''j^^2a-2'?+2 2S',(p,-2g + 2), 

' (V.9) 

Cr2,,s2p+i'(0) = 7r''^^2a-2''2S'^(p,-2g). 

The absence of terms with primeless is a consequence of the location of the real 
zeros of (r. 

2. OddD-M = 2q + l 

(a) M = 2p + 2 

Cr29+1xS2p+2(1) = 0, (V 10) 

Cr2'J+1xS2p+2'(0) = 0. 

This vanishing follows from known properties of Ch{x, 1/2) and, as a result of ( |V.2| ), 



W 



(b) M = 2p + 1 



1 / 1\ 2a~^'^+^ 
Cr2,+1xs2p+i(1) = vr<?+2r + - — — -E^(p,-2g + l), 

V 2/ 2p ! .^.oN 

Cr2,+1xs2p+i'(0) = TT'^+^r -g-- — — -S^(p,-2g-l). 



2; (2p)! 

When studying the four kinds of sums ( |V.4| ) and ([V.7|) for p — > 00 and finite 22; of the type 



—2q— 1, — 2g, — 2g + l, — 2g + 2, one has to consider the behaviours of the and l3k coefficients, 
which vary within the ranges 

ao(p-l) = 1, ••• ,ttp(p-l) = - 



vr ' (V.13) 
/3o(p-2) = 1, ... ,/3p-i(p-2) = [ip-m\ 

and also satisfy 



- 1) ~ - 2) ~ as p ^ 00. (V.14) 



For the Pk coefficients, this property was already observed in [[Tl|] . 
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We must also take into account the following asymptotics: 

CH(-(2n + l),l/2) C^(_(2n + 1)), 

Cj,{-{2n + l)) CR(-(2n + l))logn, 

G(-(2n+l),l/2) ~ -cy-(2n + l)), 

which are valid for n ^ 1, and follow from known results about the gamma, Riemann and 
Hurwitz zeta functions. 

We show, as an example, the case of —2q + 1). Including the prefactor l/{2p)\, we 

p-i 

denote the terms in that sum by S/3(p, — 2g + 1) = (— 1)^^"^ ^ t^, where, as one may check for 

k=0 

large p, < to < • • • < tp-i. Therefore 

|S^(p,-2g + l)| <p Vi, (V.16) 

Combining the information we have with the Stirling approximation for the factorial (or F) 
functions, we get 

2(2q^ [(p-l)f (2g+l)!C(2g + 2) p-V^ 
P Vi = P (2^)2,+2, C(2g + 2)-^^^yi 2^^+%^^+3/2 

As a result, 

S^(p, — 2g +1) 0, for any finite positive q, 

(V.18) 

p ^ oo. 

The other three types of sum have the same property but, since the proof is of similar nature, 
the details are omitted. In consequence, 

Cro-Mxtm(I) 0, 

'(0) ^ 0, (V.19) 
M oo, finite {D - M), 

i.e. the two-point Green function at equal points and the Tr log contribution to the free energy 
vanish in this case of the D —>■ oo limit. The importance of that limit lies in the chance of using 
the spacetime dimension as a perturbation parameter in field theory, with the advantage that 
it is then possible to obtain nonperturbative results in the couphng constants , such as Green 
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functions for quantum fields in the Ising limit. In this spirit, expansions in inverse powers of 
the dimension have proven quite useful in atomic physics 

The m > case is mathematically more involved. A possible way out is the construction 
of a power series in am by combining the preceding results with a simple binomial expansion. 
Such a method leads to 



Cro-«xs"(s) = 2a 



2s~D+M 



TT 2 



E 



am?''r is + k 



T{s) k\ 



D-M 



X 



S„(p,2s + 2A;-L) + M), forM = 2p + 2, 
S^(p,2s + 2A;-L) + M), forM = 2p+l. 



(V.20) 



V-A A calculation in x S' 



For the degeneracy of each spherical mode (see ([A.14|) ) is just (/ + 1)^. After applying 
standard Mellin-transform techniques, we end up by writing the zeta function as follows: 

^3/2^2{s-l) 



Cr2xS3(s) 



Tis] 



--(am' 
2^ ' 



f s - ^) - ^4°^ (s - 1, am) + {s - 1, am) 



(V.2i: 



where 



(V.22) 



d 

'MV^5"; ^ '^^ rft^ (t) ■ 

Although not exactly like (|1V.3|) , these integrals may also be written as Dirichlet series involving 
modified Bessel functions. Furthermore, for M = 3 such Bessel functions are expressible by 
finite series, and it is then possible to interchange the summations finally arriving at finite sums 
of polylogarithm functions. In this way, we get 



Cr2xS3(1) 



--[am 



,3/2 



Liofe" 



-2-Kam\ 



[am] 



+7r2amLi_2(e"'™") 



Cr2xS3'(0) =A'k^ 



[am 



,5/2 



[am] 



\ (27ram)^ 
'Li3(e~2^'^™) 



+ 



15^ ' ^ ' V (2'/ram)^ 
+7r2(am)='Li_i(e-2™™)] . 



+ 3 



Li3(e-2™™) 
(27ram)^ 

m('g-2™m> 



Li5(e 



-2ivam\ 



{2T[am)^ 



{2'nam) 



(V.23) 
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Since we are in D = 5, Ac must be the value of A satisfying the second gap equation for 
6 = 0, m = 0, i.e. 



1 



— = G{x,x;0,g) 

Ar 



1 



(27r)2a' 



m=0 



Li3(l) 



1 



^(3). 



Next, we replace ^ with this critical value and solve numerically the second gap equation for 
6 = only. The critical value obtained for am is 



(am)c = 2.2689. 

Calculating Cr2xs»'(0) at am = {am)c, we find the finite contribution to the free energy density 
at the critical point: 

V-B Application to T" x 

The zeta function in this space may be written 



27r 



-2s 



6' 



where we are using the notation 



ni,...,nM I I 



/ + — ] + a^m^ 



2na 



(V.24) 



(V.25) 



and d{P,l) is the degeneracy of the Ith P-dimensional spherical mode ( |A.14|) . In the m = 
case, using ([V.6|) this is put in the way 



X E E(^ 

ni,...,nM I 



P-1 



p-i-2fc r 



-2s A:max(P) 

5- E (-i)'A(P) 

A:=0 

-(nj + ... + nl) + 6[l + ^- 



21 -s 



(V.26) 



For odd P, is an integer and simple properties under dual transformations 5 ^ | in the 
sense of sect. 4 in the second ref. of IQ may appear, after making the replacements explained 
there. It is quite clear that it P es even, there can be no such invariance — in whatever sense. 
However, it can be achieved considering antiperiodic (instead of periodic) field solutions, which 
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lead to the change — nj + i, and making some further alterations. In particular, for M 
P = 2, rf(2, /) = 2{l + 1/2) and / + (P - l)/2 = / + 1/2. As a result, one can then put 



Ct1xs2(s) oc ^ + i 

n,l 



J "+2) +n'+2 



Then, the ensuing free energy will be invariant if one cares to replace the ordinary 5-derivative 
with the 'fractional-derivative' from [Q. In a similarl way, by adequately modifying the defini- 
tion of the free energy we may get the 5 \ invariance in higher dimensions. 
Combining the general forms of the previous calculations we find 

M-2s 



-2s 



^M/2 (22L\ 00 
1=0 



Tis) 



Im s 



27r 



m(P,/)) + (|^m(P,/)y 



T s 



M 



where 



m(P, /) 



+ IV?. 



(V.27) 



(V.2^ 



After expanding Im into a Dirichlet series of modified Bessel functions, we realize that for 
P ^ 00 this part may be neglected because it is exponentially vanishing. Calling Ct^xs^'^'^^ 
the rest, one gets 



2a2^ (f r (5 - f ) ^ j 2. - M), for P = 2p + 2, 

S^(p,2s-M), forP = 2p+l, 



Tis) 



(V.29) 



(47^)Af/2 

which is, up to a constant, the zeta function for R*^ x with m = 0, whose properties have 
already been considered. 



VI Higher-derivative 0(N) non-linear sigma model in 



It is interesting to observe that the model ( |II.2| ) may be easily generalized to have higher- 
derivative terms. In order not to have to study higher-derivative conformally invariant opera- 
tors, we limit ourselves to R'^^*-'^ x T*"^, which is relatively simple due to its conformal flatness. 
Then we may write 



Z[9] 



(vi.i: 
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(VI.2) 



Repeating all the steps in section 1, we get 

N 

The Green function and the Tr log part of the above expression are directly linked to the 

associated zeta function, which admits the following power expansion in po"^/^: 



Tr log(n^ + ^) ~ / ^^^V9 



m 



TT^-^ o-M 



02^(_l)fer(25 + 2A; 



r(.) 

D-M 



-a 4 



D-M 



D-M „ / 

2^ r( 



D-M+2 ' 
4 



22fcr (s + A; + i) 



2s + 2k- 



D-M\ f pa 



4s+Ak-D+M' 



2n 



(VI.3) 

where is the usual Epstein zeta function. This function has poles at s = ^ — A;, A; = 
0, 1,2,.. ., with the exception (if they coincide) of s = 0, —1, —2, . . ., where it is finite. 

Mathematically speaking, this object is harder to deal with than an 'ordinary' zeta function 
( i.e. one for a second-order operator). It can also be expressed by a series of hypergeometric 
functions: 



S-RD-M^rj.Myo) 



— 



2^r(s) VP/ 



r(. + i) 



X E 2i^l I « 



4 ' 4 
or by the integral representations 



D-M D-M 1 1 I pa 
S : ^ "S + 7:; - 



1/4- 



(V1.4) 



2n 



in') 



2\-2 



[n 



-2s 



, . TT 2 era ' 4 



X 



Cr,£>-m''x'j,m(s) — 

rl/4\ 



271 



D-M-Gs-l 



1d,m I s. 



per 



27r 



+ 



pa 



23.-l/2r(s) 

1/4 \ f-es-i 



r(4. + i-f) 



27r 



2^-3«r (-2s + 2 + f) 



(VI.5) 



where 



1=1 



fC 



dt t''-^Js+i/2{aH) e-Tl"l 



V ; ; 

ne(N*) 

Here J is the first species Bcsscl function. 

A difi^erent strategy is to regard + cr as (□ + icr^/^)(n — ia^^"^). Then, 



(VI.6) 



1 1 r 

-Tr log(n Vfj) = :t Tr log(n + i(7^/^) + Tr log(n - ia^^^) 



P_ 
27r 



D-M 



.(□+i(TV2) 



Re C^^.^ (0). 

(VI.7) 
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Taking advantage of the calculation in sect. for — □ + when D = 2d + 3, we arrive at 
the following expression for the finite part of ^(pa^/^): 

( d+l 



1 



(47r)'i+3/2 



} nd+3 /-Sr {d+l + k)\ , i/4,rf_|.i_fc 



r d+l-k 

i 2 Li 



jl/2p^l/4 



(VI.8) 

Observing that the values at the critical point have to coincide with the extremals of W/N as 
a function of pcr^/^, these points are found by numerical examination of ( |V1.8| ). 



d 


min of -2 ) CR^-^+axTi (0) 




1 


1.65 


-0.1459 


2 


2.36 


-0.1088 



For d =0,3 and 4 no solution has appeared, while for c/ = 5 we get 3.03 and -0.3141. Comparing 
with numerical estimates for the model ( [11. 2| ) on the same background, we see that the properties 
of the higher-derivative model (|VI.1| ) are drastically different. 



VII Conclusions 

In the present paper, using zeta-function regularization, we have calculated the free vacuum 
energy (or Casimir energy) for the D-dimensional 0{N) non-linear sigma model on some spaces 
with constant curvature. Explict expressions for the free energy have been obtained at the 
critical point (when the gap equations have been used) and also in non-critical regime (when the 
gap equations have not been used). For all those spaces with D = 3,4, the explict expressions 
for the free energy may be easily used in studies of the quantum antiferromagnet |18[. Moreover, 



for 2 < D < 4 it is known that the 0(iV)-model in R-^ possesses an order-disorder phase 



transition and hence its free energy may be very useful for studying dual properties. In this 
respect, the generalization to curved backgrounds is of interest. 

In particular, we have shown (see eq. ( |V.10| )) that the free energy of the model in R^^"*"^ x 
g2p+2 vanishes. This generalizes the corresponding result of ref. |P in R^ x S^, and shows that 
in D dimensions the free energy at the critical point for R^''"'"^ x S^^^^ is the same as that for 
R^, in accordance with the conformal equivalence between both manifolds. 
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From another viewpoint, for the development of connections of two-dimensional conformal 
field theory (which is a very powerful tool) with higher dimensional models, it is important 



to study the modular properties of such theories 0]. In subsect [V-B| we have checked dual 
invariance in Cardy's sense for m = and antiperiodic modes in x S^. Hence, the higher- 
dimensional 0(A^)-model may serve as a very useful toy model for studying dual symmetries, 
which have become quite popular in recent studies of strings and supersymmetric QCD. 

Finally let us note that the D-dimensional 0{N) non-linear sigma model may be consid- 
ered as a toy model of quantum field theory in a Kaluza-Klein framework for studying the 
question of spontaneous compactification, for example. As usual, Kaluza-Klein theories are not 
renormalizable. However, in the model under discussion we may still use the -^-expansion in 
order to control somehow the quantum corrections. 

Observe also that in addition to the -^-expansion one can calculate the free energy as an 
expansion in inverse powers of the dimension. Such calculations in frames of the toy model 
under consideration may be useful in atomic physics M. 



A Appendix: on the calculation of zeta functions 
I-A General considerations 

Let M he a spacetime of the type R^"^^ x , D > M, with the operator —O + m'^ + ^R 
acting on the whole manifold. Since — □ = —O-^d-m — Oj^m^ the spectrun has the form 
+ Xn + + ^R, p G R-^"^"^, A,„ e Sp (-□_yv,jv). Then, the global zeta function on = 
^D-M ^ defined as 

CRi.-M,^M(s) = / d^-^p J2 ip' + An + + ei?)-^ (A.l) 

n: 

A„gSp(-D^j,,^) 

Such a definition is purely mathematical, i.e. the usual physical factors coming from state 
normalization are here absent. Notice that each momentum integration is introducing a further 
mass dimension. The arbitrary mass scale fi, typically supplied in order to render the function 
dimensionless, is not present either. All these elements may be included at a later stage. 

After doing the p-integrations, this zeta function can be written in terms of the one for 
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, , £^r(s-^) / D-M\ ,^ , 

C,^o^m^^m{s) = TT 2 ^Ca^a^ (^S ^ j , (A.2) 

with 

Ca,m(^)= ^ (A„ + m2 + ei?)-^ (A.3) 

n: 

where possible degeneracies must be accounted for into this sum. 

First, we make the following hypothesis: only the T function has singularities at s — 
s = 0, 1, — and, obviously, these poles can be encountered for even D — M only — while Ca^aj 
and its derivative are finite at these points. Such an assumption is right when M.^^ = T^^ 
for odd or A^^^ = S*^ for the massles case. Using the expansion of the F functions, both 
around their poles and around regular points, we obtain the generic results: 



1. Even D-M = 2q 

1 



Cr29xa^a^(1 + e) = vr^ 



(g-1)! 

(-1)9 

Cr29xa^m'(0) = ^"^^ {[^(^) + ^] i-q) + (mm (-g)} • 



2. OddL>-M = 2g + l 



Ca,m {-q + 1) + C^m' {-q + 1)| + 0(e), 



(A.4) 



Cr29+1xA4a^'(0) = 7r«+^F (^-q - Ca4" (-9 - ^ 



(A.5) 



When, on the contrary, Cm'^'{^) has poles at 2; = s — s = 0, 1, this function has to be 

Laurent-expanded around these singularities before taking the limits s — 0, 1, and a somewhat 
different calculation is required in every particular case. 



I-B Zeta function on 

We study the torus only. Once we have the zeta function on T*^, we can find the one on 
Cn^-My^j^M by application of ([A.2|) . So, when looking at Ct^n{z) we will have in mind the 
arguments z = s — ■, s = 0, 1. 
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For simplicity we take all the radii equal and with value p. Then 



27r 



-22 



Ctm{z) = l—j Zm [z, 



pm\ 
27/ 



(A.6) 



where 

oo 

ZM{z,a)= {nl + ... + nlj + a^)-' (A.7) 

ni,...,nj(/=— oo 

is an Epstein zeta function with inhomogeneous term. After Mellin-transforming we obtain, 

-2z 

I -ZTT \ 

Ctm[z) 



27r 
. P , 



Iz) Jo 



Viz) 



The Jacobi theta function 6{x) — given by ( [IV. 4|) — has the property 



9{x) 



6 

X \x 



(A.9) 



Taking advantage of this, we separate the part which diverges at t = — i.e. the n = 
contribution in the 9 function — and write 

/ 

/o 



27r\ 




P J 




27r\ 




P ) 


r(.) 



TT 



1 



V27/ 



M~2z / M 

T Z- — 



M~2z 



r z- 



M 



(A.IO) 

Now the t-integral, that (in accordance with ( p.V.3| )) we have called Im (-^^j^)) contains no 
small-t singularity. This expression has the additional advantage that its second term exhibits 
the poles of this function at z = ^^Y^^ ■ ■ ■ After writing 9^'^ as a binomial expansion (in the 
n = term of 9 and the rest of the summatory), we use the integral representation 



/ dt f'-'e 
Jo 

and end up with the Dirichlet series 



u/2 



K^{2Vab) 



{A.ll] 



M-2z J£ / 

I 



These results, combined with ( [A.2| ), give rise to ( [iV.2| ),( [rv.5D for CR-f'-"xT"(s). 



n 



^+1 g KM/2-z{pm^nl + .. 

ni,...,n,=i {pmJnf + ... + ra,2)AV2-2 



(A.12) 
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I-C Zeta function on S 



M 



The same method is employed to derive the expressions in Cn^-^'^M^i^) from the ones for S^^. 
Taking into account the known spectrum of the D'Alembertian on S^^, the Riemann curvature 
on this space and the conformal couphng value of ^, we readily obtain the eigenvalues of our 
operator and construct its zeta function 



CsM{z)=Y.d{M,l) 



1=0 



M- 1 



+ m 



with degeneracies (see e.g. p^ ) 



(A.13) 



l + M-2 \ 1 + ^ 



(A.14) 



Setting m = one obtains 

,2z oo / 



2a 



M 



1=0 



2a 



2z P 



l + M-2 \ ( M-l\^-2^ 

/ H 

/ IV 2 



(2p + 



— 5](-l)'=afc(p-l)CH(^-l + 2^-2p + 2A;,p + -j, for M = 2p + 2 



A;=0 

-— 5:(-1)'=A(p-2)Ch(2z-2p + 2A;,p) 



for M = 2p + 1, 



(A.15) 



with the Ofc and coefficients as written in (|V.5|) . The Hurwitz functions in ( |A.15|) may be 
reexpressed with the help of the identities 



C«(^.i)-|:(« + i)"'.C«(4) = (2-i)CBW. 

n=\ 



(A.16) 



Afterwards, taking advantage of the properties 

E(-i)'«'^(p-i)(^+2) 

E(-im(p-2) (n + if (^-^-'^^ 

k=0 



0, when < n < p — 1, 
0, when < n < p — 2, 



(A.17) 
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we may put 



2a^^S„(p,2z), for M = 2p + 2, 
2a2^E/3(p, 2z), for M = 2p + 1, 



(A.18) 



where the Ea(p, 2^;) and 2^;) are the ones defined by ( |V.4| ). Observing those finite sums, 

one locates the singularities of (s^i{z), which come from the pole of (h{x, a) or (r{x) at x = 1: 



Cs^i{z) is finite and we may therefore apply ( |A.4|) , ( |A.5| ), thus getting the values of Cr-d-a/xs*i^(1)) 
CRO-MxsAf'(O) given in sect. 0. 
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